Quantum shutter approach to tunneling time scales with wave packets 
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The quantum shutter approach to tunneling time scales (G. Garcfa-Calderon and A. Rubio, Phys. 
Rev. A 55, 3361 (1997)), which uses a cutoff plane wave as the initial condition, is extended in 
such a way that a certain type of wave packet can be used as the initial condition. An analytical 
expression for the time evolved wave function is derived. The time-domain resonance, the peaked 
structure of the probability density (as the function of time) at the exit of the barrier, originally 
found with the cutoff plane wave initial condition, is studied with the wave packet initial conditions. 
It is found that the time-domain resonance is not very sensitive to the width of the packet when the 
transmission process is in the tunneling regime. 

PACS numbers: 03.65.Xp, 03.65.Ca, 03.65.Ta 



I. INTRODUCTION 



Tunneling is one of the most important quantum phe- 
nomena that has been widely applied in science and tech- 
nology. For years, the stationary treatments of tunnel- 
ing were sufficient for many practical purposes, and the 
details of tunneling dynamics were not urgent issues to 
investigate. This is not the case anymore. The inter- 
est in tunneling dynamics is increasing as, for example, 
the number of carriers involved in tunneling events de- 
creases due to the rapid downsizing of semiconductor de- 
vices. In principle, the tunneling dynamics can be com- 
pletely understood if one can solve the time dependent 
Schrodinger equation, taking other degrees of freedom 
into account that affect the tunneling particles. This is, 
however, a difficult task in general. It is thus important 
to study the time scales of tunneling dynamics (tunneling 
times) in simplified models and use them for qualitative 
understanding of the tunneling dynamics in realistic sys- 
tems. There are many approaches to define or measure 
the tunneling times in simplified models Q, 0, 0, 3 , and 
the quantum shutter approach [5| is one of them. The 
present paper concerns a generalization of the quantum 
shutter approach. 

Let us consider a one dimensional scattering problem 
where a particle is incident on a potential V(x). The 
time scales that characterize the tunneling dynamics are 
"embedded" in the wave function. To extract them from 
the wave function, the quantum shutter approach |5| uses 
a cutoff plane wave as the initial condition and monitors 
how the probability density changes in time at a spec- 
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ified position (e.g., at the exit of the barrier) or in a 
spatial region (e.g., in the well region in a double barrier 
structure) to find out the time scales that characterize 
the transient behavior of the wave function (from non- 
stationary to stationary) . The use of a cutoff plane wave 
initial condition can be understood as an analogy to the 
use of a step input in the study of the temporal response 
of an electrical (e.g., RCL) circuit. 

Studies of tunneling dynamics involving cutoff plane 
waves go back back to Stevens Q , who argued the signal 
velocity under the barrier by applying the contour defor- 
mation technique developed by Brillouin jjj. This tech- 
nique allows one to decompose a wave propagating in a 
dispersive medium into three parts: the fore-runners, the 
monochromatic part oscillating with the same frequency 
as the source, and the after-runners; the signal velocity 
is then defined as the velocity with which the monochro- 
matic front moves. The under-the-barrier signal veloc- 
ity found by Stevens was, however, questioned later by 
Teranishi et al 0, Jauho and Jonson and Ranfagni 
et al jl(| ; these authors showed that the monochromatic 
front to which the signal velocity is attributed is in fact 
not appreciable in magnitude. Through these works, it 
was recognized that it is not appropriate to focus only 
on the monochromatic part of the wave. Biittiker and 
Thomas an d Muga and Biittiker thus studied not 
only the monochromatic part but also the fore-runners in 
detail, whereas Brouard and Muga turned their at- 
tention to the total wave function (under a cutoff plane 
wave initial condition) and studied its properties with an 
exact analytical expression for the wave function which 
they derived. Independently, Garci'a-Calderon and Rubio 
derived an exact analytical expression for the wave 
function with a cutoff plane wave initial condition and 
applied it to the analysis of the transient behaviors in 
the tunneling dynamics. It may be said that these stud- 
ies formed a new area of research in the field of quantum 
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dynamics, where one explores the tunneling dynamics, 
especially its transient behaviors, by using exact analyt- 
ical expressions for the wave functions with cutoff plane 
wave initial conditions. The quantum shutter approach 
and the approach by Brouard and Muga 01 use the 
same tool called M function [see Eq. to express the 
wave functions in analytical manners. They differ, how- 
ever, in the following respect: in addition to the M func- 
tion, the quantum shutter approach uses the resonant 
eigenf unctions, while the approach by Brouard and Muga 
uses an entire function that arises from a pole expansion 
of a function involved in the momentum eigenfunction 
expansion of the wave functions. In |l4j, an analytical 
expression for the transmitted wave was obtained in sole 
terms of the M function and the poles and the residues 
of the transmission amplitude. Another interesting and 
unexpected aspect of the quantum shutter approach is 
that it has a close relationship with the consistent his- 
tory approach to the tunneling time problem as shown 
in 14] . In particular, the probability density at the exit 
of a rectangular barrier under a cutoff plane wave initial 
condition, a quantity of major concern in the quantum 
shutter approach, coincides, when properly normalized, 
with a function G p (t) that is defined in the consistent his- 
tory approach. Function G p (t) allows us to associate the 
transient behaviors of the wave function with the interfer- 
ence between Feynman histories with different tunneling 
times, which provides a novel viewpoint to the transient 
behavior. In this way, the quantum shutter approach is 
related not only to the stream of research that began 
with the work of Stevens but also to a relatively new ap- 
proach, the consistent history approach, to the tunneling 
time problem. In Ref. [lq. the consistent history ap- 
proach was used to give different tunneling times in a 
unified manner. 

The quantum shutter approach, formulated for a gen- 
eral (but finite-ran ge) potential, has been applied to sin- 
gle barriers 0,0,^3; double barriers and to a more 
general superlattice structure 0] . For the case of tunnel- 
ing through a single barrier, the probability density at the 
exit of the barrier was found to have a peaked structure, 
which is called the time-domain resonance. The time i p , 
corresponding to the peaked value of the resonance, has 
been studied in detail as a function of the system param- 
eters and has been proposed as a tunneling time scale 
|18|. For the case of tunneling through a multi-barrier 
resonant structure, the buildup process in the well re- 
gion^) has been also studied, with the expectation that 
it has important implications on the speed of resonant 
tunneling devices [{J [Mm- These studies, while demon- 
strating the usefulness of the quantum shutter approach, 
raise a natural question: To what extent do the results 
obtained from the approach, in particular the existence of 
the time-domain resonance, depend on the special form 
of the initial conditions? To answer this question, it is 
necessary to extend the approach to a more general class 
of initial conditions. 

The purpose of the present paper is to provide an ex- 



tension of the quantum shutter approach to a type of 
wave packet initial condition, which reduces itself to the 
cutoff plane wave initial condition in the limit of an infi- 
nite packet width. This extension does not change the an- 
alytical scheme of the approach, so that the time evolved 
wave function continues to be available in an analyti- 
cal form. We apply the extended approach to tunneling 
through a square barrier to study how the time-domain 
resonance is affected by the width of the initial wave 
packet. 



II. FORMALISM 

The quantum shutter approach assumes an arbitrary 
but finite-range potential V{x) [in this paper V{x) is such 
that it vanishes for x < and for x > d] and uses, for 
example, the following form of the cut off plane wave as 
the initial wave function: 



#(x,0) = 



2i sin fco^ 




for x < 
for x > 0, 



(1) 



where ko is the wave number. This setup corresponds to 
the physical situation where a beam of particles with en- 
ergy Eq = fi, 2 fcg/2m (though not exactly monochromatic 
due to the sharp front of the wave) impinges on a shutter 
placed at x = 0, just at the left edge of the potential; the 
tunneling process begins with the instantaneous opening 
of the shutter at t — 0, enabling the incoming wave to 
interact with the potential for t > 0. When the initial 
condition is given by Eq. Q), the exact solution of the 
time dependent Schrodinger equation along the transmit- 
ted region x > d is found to be [Til fl6l| 

(a;, t) = T(k )M (x, ko\t)- T(-k )M (x, -k ; t) 



-2fco J2 



n— — oo 



U2 _ h,2 



M(x,k n ;t), (2) 



where k n is the n-th pole of the transmission amplitude 
T (poles lie in the lower- half of the complex k plane), r n 
is the associated residue, and the M functions are defined 
by 



M(x,q;t) = — 



dk- 



a ikx—ihk t/2m 



; e ( '" 2/2t ^(i!/,), 



(3) 



(4) 



where q = fe„,±fco; the function w(z), which is often 
called Faddeeva function, is related to the complex com- 
plementary error function as w(z) = e~ z erfc(-iz) |20| . 
and y q is given by 



y q = e-™ /4 



hq 



(5) 



The w function appears in many fields of physics and 
mathematics, so that it has been well studied and its 
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properties have been well understood |2l| . Some com- 
puter programs are available for numerical calculation of 
the w function 

Let us now consider the following initial condition: 



*(x,0) = A 



dk 



ik.r 



k — &o + iA 



(6) 



where A = y/A{l + (A//c ) 2 }/27r with A > 0, and "c.c." 
stands for complex conjugate. An important feature of 
this 0) is that it automatically vanishes for x > 0. 
This is immediately seen from the fact that the integrand 
e lkx I (k — ko + iA), which corresponds to the Lorentzian 
momentum distribution centered at hko with width HA, 
has a simple pole only in the lower-half of the complex 
fc-plane. An explicit expression for ^(x,0) can be easily 
obtained by the method of residues. We find 



*(*,0) = {^ 



■ Ax sin h 



Q.V 



for x < 
for x > 0. 



(7) 



This represents a wave packet. A measure of the packet 
width is 1/ A. We can easily prove that the wave packet is 
normalized, i.e., J dx\*f?(x, 0)| 2 = 1. If the above ^(x,0) 
is multiplied by a constant i/V~A and the limit A — * is 
taken, Eq. (JTJ is reproduced. The above ^>(x,0) is there- 
fore a wave packet counterpart of the cut-off plane wave 
initial condition. Equation © thus leads us to a natu- 
ral setup of the shutter problem with a normalized wave 
packet initial condition which vanishes automatically for 
x > due to the Lorentzian momentum distributions. 

Let us derive an analytical expression for ^>(x,t) under 
the initial condition given by Eq. |7jl. For definiteness, 
we shall limit our attention to the analytical expression 
only in the transmitted region x > d. We start from the 
following expression for the time evolved wave function 
for the transmitted region: 



dk 



<p(k)T(k)e 



ikx — ihk t/2m 



(8) 



where 4>{k) is the fc-space wave function (i.e., the Fourier 
transform of the initial wave function) defined by 



dx 
■ ( 



-ikx 



*(x,0). 



(9) 



Equation JSJ follows directly from the eigenfunction ex- 
pansion of the wave function in the transmitted region. 
Substituting Eq. © into Eq.©, we have 



(j>{k) = V2ttA 



1 



1 



k - k + iA k + ko + iA 



(10) 



Next, we expand the transmission amplitude in terms 
of its complex poles and the corresponding residues by 
using a special form of the Mittag-Leffler theorem due to 
Cauchy |2^] . It may be expanded as 0, |24| 



T{k) = 



OO 

E 

n— — oo 



(ii) 



The substitution of Eqs. (|TU|l and (fTTft into the right- 
hand side of Eq. |JHJ yields the following quantity: 



1 



1 



k - k + iA k + ko+iA 



1 1 

k kn. km. 



which, upon expansion, gives four terms. To two of the 
four terms, we apply the partial fraction expansion 



1 



1 



1 



k ± ko + iA k — k„ 
1 



±fco + k n + iA 
1 



k — k n k + ko + iA 



(12) 



and express the resultant k integrals in terms of M func- 
tions. We then have 



x E 

n 

-E 

n 

-E 



-2niA 



't - K n 



iA 



k n + iA ko — k n — iA 



M(x,-k ~ iA;t) 
M(x,k n ;t) 



(13) 



Due to Eq. (|llf> . the first and the second sums over n 
in the square brackets in Eq. (|13fl give T(ko — iA) and 
T(~ko — iA), respectively. We thus arrive at 



*(z,t) 



i^A{l + (A/k ) 2 } 
T(ko-iA)M(x,k -iA;t) 

- T(-ko - iA)M(x, -k - iA; t) 



2k J2 



K - (k n + iA)' 



:M(x,k n ;t) 



(14) 



This is the analytical expression for the time evolved 
wave function in the transmitted region under the ini- 
tial condition given by Eq. (JJJ. The analytical solution 
under the cutoff plane wave initial condition, Eq. can 
be correctly reproduced from Eq. (|14|) if we multiply Eq. 
(|14fl by a constant i/y/A and then take the limit A — > 0. 
It is also possible to derive an analytical expression for 
^(x, t) in terms of M functions in other regions of space, 
although we concentrate on the transmitted region. 

To use Eq. (|14|) . we have to find the poles {k n } and 
calculate the residues {r n } and the w functions numeri- 
cally. The residues may be calculated in general by using 
the simple relationship [l4| . 



r n = iu n (0)u n (d)e ik " d , 



(15) 
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where, as discussed in the appendix A of Ref. the 
resonant eigenfunctions u n (x) are solutions to the time- 
independent Schrodinger equation 



d 2 u n (x) 
dx 2 



i 2 2rn 
k n - —^-V{x) 



Unix) = (16) 



satisfying the outgoing boundary conditions, 



dx 



u n (x) 



-ik n u n (0); 



dx 



u n (x) 



and the normalization condition, 

u 2 n {0) + u 2 n {d) 



u n {x)dx + i- 



U n (d), 

(17) 
(18) 



For a rectangular potential of height Vq and width d, the 
resonant eigenfunctions read, 



u n (x) =C n [e i *» x +D n 



{0<x<d) (19) 



where q n = [k% - k^} 1 / 2 , k v = y/2mV /U, D n = 
(q n + k n )/(q n — k n ) and C n may be obtained from the 
normalization condition given above. Alternatively, one 
may use the following explicit relationship between k n 
and r„ to calculate r„ : 



k 2 v fl 2 e 



ik,, d 



ky {krid - 



2i) sm(dy/k% 



(20) 



One can derive Eq. (|20() directly from r„ = 
exp(—ik n d)/g'(k n ), where the prime stands for d/dk and 
g(k) = exp(— ikd)/T(k); the exact analytical expression 
for the transmission amplitude T(k) is available in the 
standard textbooks of quantum mechanics. 

Readers might have noticed that both Eqs. and Q 
correspond to an initially vanishing probability current 
density, i.e., J(x,0) = 0. This is, however, not a general 
feature of the quantum shutter approach. The approach 
can be formulated even if $(x, 0) 



2ismk x = e %k ° x 



-iknx 



in Eq. Q is replaced by \&(x, 0) = ae lk ° x + be 



-ik\)X 



with arbitrary constants a and b. For this general plane 
wave initial condition, J(x. 0) 7^ in general. In the same 
vein as above, we can construct a wave packet counter- 
part, for which J(x,0) 7^ as well. For this wave packet 
initial condition with nonzero flux, one may also obtain 
an analytical expression for the time evolved wave func- 
tion in terms of M functions. 



III. EXAMPLE 

In the rest of this paper, we apply Eq. i|14|) to tun- 
neling through a rectangular barrier to study how the 
time-domain resonance [lj| depends on the width of 
the incident packets. Assuming a rectangular barrier 
of height Vq — 0.3 eV that extends from x — to 
x = 4 nm, we have calculated the probability density 
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FIG. 1: Plot of |*(d, t)\ 2 at the barrier edge x = d = 4 nm as 
the function of time in units of the free passage time t{. The 
initial condition is given by Eq. Q in (a) and by Eq. JZj 
in (b) and (c). We calculated numerically 1000 poles in the 
3rd quadrant and also in the 4th quadrant in the complex k 
plane, and used them with Eqs. and 1141 to plot these 
graphs. 



at the exit of the barrier, \^(d, t)\ 2 , for a particle of ef- 
fective mass m = 0.067m o (m c being the bare electron 
mass) with (central) energy Eq = h 2 k 2 l /2m = 0.01 eV 
(fco ~ 0.133 nm -1 ). The above refers to typical semicon- 
ductor heterostructure parameters 0|. Figure H shows 
the results for different values of A, where the time axis 
is in units of the free passage time U = md/hko « 17.5 
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fs. Figure [I] (a) shows the plot of \^(d, t)\ 2 for the case 
of A — calculated from Eq. J2J); we can see a tran- 
sient behavior of the probability density that starts from 
zero, increases to the maximum, and then approaches its 
asymptotic value |T(fco)| 2 as shown in the inset of Fig. 
1 (a). The peaked structure is an example of the time- 
domain resonance |16|. In the present example, the time 
t p that gives the maximum of \^(d, t)\ 2 is about 0.303tf. 
Figures (b) and (c) obtained from Eq. (jT4*|) show re- 
spectively the case of A/fc = 0.75 (A w 0.10 ran -1 ) and 
the case of A/fc = 4.0 (A ss 0.53 ran -1 ). Unlike the case 
of Fig. [I] (a), the probability density approaches zero for 
large values of t in Figs. ^ (b) and (c). This is simply 
because the particle is eventually reflected or transmitted 
for the wave packet initial condition. 

We first find from the graphs that the peaked nature 
of \^(d, t)\ 2 is not lost (i.e., the time-domain resonance is 
preserved) even for relatively large values of A (and thus 
for relatively narrow packets in configuration space). In 
the case of Fig. 12(b), the packet width 1/A is more than 
twice the barrier width. In the case of Fig.^(c), where 
the time domain resonance is still clear, the packet width 
is less than half of the barrier width. 

Secondly, we find that t p does not depend significantly 
on the value of A in the tunneling regime. In the case 
of Fig. ^ (b), where the average energy of the particle, 
(E) = {1 + (A/k ) 2 }E , is approximately 0.016 eV, we 
have tp/t{ » 0.288, so that the shift of t p /t{ from the case 
of Fig. (a) is only about 5%. In the case of Fig. ^ ( c ), 
where (E) « 0.17 eV, we have t p /t{ « 0.221, which shows 
as much as 27% shift from the case of Fig.n(a) . Although 
the average energy is below the barrier height in both 
cases, we must note that only the case of Fig. Q] (b) can 
be associated with tunneling as explained below. 

In what follows we refer to two different approaches 
to distinguish between tunneling and non-tunneling pro- 
cesses. The first one involves a stationary analysis 
whereas the second one deals with a time-dependent de- 
scription. For wave packet initial conditions, the first 
approach relies on the computation of both the tunnel- 
ing probability P un der (i-e., the probability of under-the- 
barrier transmission), and the non-tunneling probability 
■fovcr (i-e., the probability of over-the-barrier transmis- 
sion) defined by 



under — 



dk\T(k)\ 2 \cf)(k)\ 2 



P OVBr = / dk\T{k)\ 2 \^{k)Y 

>k v 



(21) 



(22) 



In the case of Fig. Q (b), -Punder ~ 0.00161 and P OVC r ~ 
0.00117, so we affirm that under-the-barrier transmission 
slightly dominates, and in this sense, the transmission 
is in the tunneling regime. In the case of Fig. ^ ( c )i 
-Punder ~ 0.0111 and P ver ~ 0.0426, so that over-the- 
barrier transmission dominates and the process is not in 
the tunneling regime. The second approach deals with 



a time-frequency analysis |25j, where we introduce the 
local average frequency o; av (the instantaneous frequency 
of the wave function) and the instantaneous bandwidth 
a (the spread of frequencies around w av ) as follows pEj : 



Wav(i) 



a{t) 



Tin I if 

* dt 



Re 



1 a* 

^~dt 



(23) 



(24) 



The local average frequency was used in Refs. 0, 
and the local bandwidth in Ref. in the studies of 
tunneling time. It is immediate from Eq. I|24|) that 
(7 = 0, i.e., the wave function has a single instantaneous 
frequency, when d\^>(x, t)\ 2 /dt = 0, which holds at the 
time domain resonance peak at t = tp. Therefore, the 
time domain resonance peak is characterized by a sin- 
gle instantaneous frequency a-> av , or equivalently, by a 
single energy 7ia; av . This was exemplified in Ref. |18| . 
We have calculated w av at the time domain resonance 
peak, namely, w av (i = t p ), to find that w aT /wy , with 
lu Vq = V /h, is ps 0.792, 0.944, and 1.583 for the case of 
Fig. 12(a), (b), and (c), respectively. This implies that 
the time domain resonance peak is associated with tun- 
neling for the case of Figs. ^ ( a ) and (b), but it is not for 
the case of Fig. [2(c) since u av /ujv > 1. We have thus 
shown, using two different approaches, that Figs. ^ (a) 
and (b) correspond to tunneling, but Fig.Q](c) does not. 

As shown above, t p is not very sensitive to the width of 
the packets in the tunneling regime. This in turn justifies 
the use of the cutoff plane wave initial condition as far as 
the estimation of the value of t p in the tunneling regime 
is concerned; t p is expected to characterize the earliest 
tunneling response of the system and thus, it would 
be of relevance for device applications. A comparison of 
t p with other tunneling time scales, such as the delay 
time, the Biittiker traversal time, and the semiclassical 
time, has been given elsewhere [lq . 

In Ref. [16|. the time-domain resonance was also 
studied with a cutoff pulse initial condition given by 
^(x,0) = 2isinfcoa; for — a < x < 0, and zero other- 
wise. It was found that the resulting time-domain res- 
onance structure is almost identical to the one obtained 
with the semi-infinite cutoff plane wave initial condition, 
i.e., Eq. (JTJ. This is consistent with our result that the 
time-domain resonance structure is not quite sensitive to 
the packet width. On the other hand, the oscillatory be- 
havior of the probability density in the "post resonance" 
time domain (see Fig. 7 in Ref. ^(|) cannot be seen 
in our results. This implies that the probability density 
in the "post resonance" time domain is sensitive to the 
shape of the incident wave. 
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IV. SUMMARY 

In summary, we have extended the type of initial con- 
ditions used in the quantum shutter approach from shut- 
tered plane waves to a certain class of wave packets. This 
makes it possible to study the "size effect" of the packets 
on the various results that had been obtained from the 
quantum shutter approach with shuttered plane wave ini- 
tial conditions. We have derived an analytical expression 
for the time evolved wave function, Eq. JT2J, under the 
wave packet initial condition given by Eq. JjJ. Focusing 
on the size effect on the time domain resonance, we ex- 
emplified that (i) the time-domain resonance structure is 
present even when the packet width 1/ A is much shorter 
than the barrier width, a situation that is quite different 
from those considered in the original quantum shutter 



approach, where the incident wave is semi infinite, and 
(ii) the time t p at which the time resonance peak occurs 
is not very sensitive to the packet width when the trans- 
mission process is in the tunneling regime. In this sense, 
the time domain resonance is robust against the change 
of the initial conditions. 



V. ACKNOWLEDGMENTS 

N.Y. thanks the Department of Physics, UNAM for 
their hospitality. He thanks H. Yamamoto for encourage- 
ment. N.Y. carried out a part of his numerical calcula- 
tions on SX7 at the Information Synergy Center, Tohoku 
University. G.G-C and J.V acknowledge partial financial 
support of DGAPA-UNAM under grant No. IN108003. 



[1] E.H. Hauge and J.A. St0vneng, Rev. Mod. Phys. 61, 917 
(1989). 

[2] R. Landauer and Th. Martin, Rev. Mod. Phys. 66, 217 
(1994). 

[3] Time in Quantum Mechanics, edited by J.G. Muga, 
R.Sala Mayato, I.L. Egusquiza (Springer- Verlag, New 
York, 2002). 

[4] V.S. Olkhovsky, E. Recami, and J. Jakiel, Phys. Rep. 

398, 133 (2004). 
[5] G. Garcfa-Calderon and A. Rubio, Phys. Rev. A 55, 3361 

(1997) . 

[6] K.W.H. Stevens, Eur. J. Phys. 1, 98 (1980);J. Phys. C 
16, 3649 (1983). 

[7] L. Brillouin, Wave Propagation and Group Velocity (Aca- 
demic, New York, 1960). 

[8] N. Teranishi, A.M. Kriman, and D.K. Ferry, Superlat- 
tices Microstruct. 3, 509 (1987). 

[9] A. P. Jauho and M. Jonson, Superlattices Microstruct. 6, 
303 (1989). 

[10] A. Ranfagni, D. Mugnai and A. Agresti, Phys. Lett. A, 

158, 161 (1991). 
[11] M. Biittiker and H. Thomas, Ann. Phys. (Leipzig) 7, 602 

(1998) ; Superlattices Microstruct. 23, 781 (1998). 
[12] J.G. Muga and M. Biittiker, 62, 023808 (2000). 

[13] S. Brouard and J.G. Muga, Phys. Rev. A. 54, 3055 
(1996). 

[14] G. Garcfa-Calderon, J. Villavicencio, and N. Yamada, 
Phys. Rev. A 67, 052106 (2003). 



[15] N. Yamada, Phys. Rev. Lett. 93, 170401 (2004). 

[16] G. Garcfa-Calderon and J. Villavicencio, Phys. Rev. A 

64, 012107 (2001). 
[17] G. Garcfa-Calderon and J. Villavicencio, F. Delgado, and 

J. G. Muga, Phys. Rev. A 66, 042119 (2002). 
[18] G. Garcfa-Calderon and J. Villavicencio, Phys. Rev. A 

68, 052107 (2003), and references therein. 
[19] R. Romo, Phys. Rev. B 66, 245311 (2002), and references 

therein. 

[20] Handbook of Mathematical Functions, edited by M. 
Abramowitz and LA. Stegun (Dover, New York, 1965), 
p. 297. 

[21] V.N. Faddeyeva and N.M. Terent'ev, Tables of values of 
the function w(z) = e _z (1 + J Q e* dt) for complex 
argument, translated from the Russian by DG. Fry and 
B.A. Hons (Pergamon, London, 1961); B.D. Fried and 
S.D. Conte, The Plasma Dispersion Function (Academic 
Press, New York, 1961). 

[22] G.P.M. Poppe and C.M.J. Wijers, ACM Transactions on 
Mathematical Software, 16, 38 (1990); 16, 47 (1990). 

[23] EG. Tichmarsh, The Theory of Functions (Oxford Uni- 
versity Press, second edition, 1939), p. 110. 

[24] N. Yamada, (unpublished) (2002); in Meeting Abstracts 
of the Physical Society of Japan, 55, Issue 2, Part2, 
23pWD-3, 227 (2000) (in Japanese). 

[25] L. Cohen, Time- Frequency Analysis (Prentice Hall, New 
Jersey, 1995). 



